Army  Research  Laboratory 


fii-  *  *=  -Tl  f  L 

r** 


Viscoelastic  Properties  of  Polymer  Systems  From  Dissipative 

Particle  Dynamics  Simulations 


by  John  K.  Brennan  and  Jan  Andzelm 


ARL-TR-4647 


November  2008 


Approved  for  public  release;  distribution  is  unlimited. 


NOTICES 

Disclaimers 

The  findings  in  this  report  are  not  to  be  construed  as  an  official  Department  of  the  Army  position  unless 
so  designated  by  other  authorized  documents. 

Citation  of  manufacturer’s  or  trade  names  does  not  constitute  an  official  endorsement  or  approval  of  the 
use  thereof 


Destroy  this  report  when  it  is  no  longer  needed.  Do  not  return  it  to  the  originator. 


Army  Research  Laboratory 

Aberdeen  Proving  Ground,  MD  21005-5069 


ARL-TR-4647 


November  2008 


Viscoelastic  Properties  of  Polymer  Systems  From  Dissipative 

Particle  Dynamics  Simulations 

John  K.  Brennan  and  Jan  Andzelm 
Weapons  and  Materials  Research  Directorate,  ARL 


Approved  for  public  release;  distribution  is  unlimited. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering 
and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to 
comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 

13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  U.S.  Army  has  an  increasing  need  for  multifunctional  polymers  with  specific  viscoelastic  properties.  Such  polymers  are 
often  present  in  nanocomposites,  lubricants,  or  biomimetic  gels  that  are  used  as  tissue  surrogate  materials.  Rational  prediction 
of  a  polymer’s  viscoelastic  properties  necessitates  the  use  of  multiscale  modeling  since  vastly  different  time  and  length  scales 
have  to  be  considered.  In  this  study,  we  will  apply  the  dissipative  particle  dynamics  (DPD)  method  that  is  a  fast  and  accurate 
alternative  to  traditional  molecular  simulations.  The  DPD  method  is  a  particle -based,  coarse-grained  mesoscale  method.  In  this 
study,  the  simulation  of  viscoelastic  properties  is  carried  out  by  applying  an  additional,  oscillatory  shearing  force.  This 
oscillatory  shear  and  corresponding  stress  allows  us  to  calculate  the  storage  (G’)  and  loss  modulus  (G”)  that  define  the 
viscoelastic  behavior  of  polymer  systems.  As  a  demonstration  exercise,  the  G’  and  G”  moduli  were  calculated  for  a  diblock 
copolymer  at  low  frequencies. _ 

15.  SUBJECT  TERMS 

viscoelastic,  polymer,  DPD,  dissipative  particle  dynamics 

17.  LIMITATION  18.  NUMBER 
OF  ABSTRACT  OF  PAGES 

UL _ 14 

standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


19a.  NAME  OF  RESPONSIBLE  PERSON 

Jan  Andzelm _ 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

410-306-4008 


16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

1 .  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

November  2008  Final 

4.  TITLE  AND  SUBTITLE 

Viscoelastic  Properties  of  Polymer  Systems  From  Dissipative  Particle  Dynamics 
Simulations 


6.  AUTHOR(S) 

John  K.  Brennan  and  Jan  Andzelm 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRD-ARL-WM-MA 
Aberdeen  Proving  Ground,  MD  21005-5069 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


3.  DATES  COVERED  (From  -  To) 

September  2007-April  2008 

5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


5d.  PROJECT  NUMBER 

AH84 

5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

ARL-TR-4647 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


11 


Contents 


Acknowledgments  iv 

1.  Introduction  1 

2.  Dissipative  Particle  Dynamics  1 

3.  Viscoelastic  Property  Calculations  4 

4.  Conclusions  5 

5.  References  6 

Distribution  List  7 


iii 


Acknowledgments 


The  authors  would  like  to  gratefully  aeknowledge  Vietor  A.  Pryamitsyn  and  Venkat  Ganesan  of 
the  Department  of  Chemical  Engineering,  University  of  Texas  at  Austin,  and  Martin  Lisal  of  the 
E.  Hala  Eaboratory  of  Thermodynamics,  Institute  of  Chemical  Process  Eundamentals,  Academy 
of  Sciences  of  the  Czech  Republic. 


1.  Introduction 


Materials  for  many  applications  relevant  to  U.S.  Army  needs,  ranging  from  lubricants  to 
nanocomposites,  are  often  designed  to  attain  required  viscoelastic  properties.  Biomimetic  gels  as 
tissue-surrogate  materials  used  in  designing  devices  to  protect  Soldiers  against  blunt-force 
trauma,  the  permselective  membranes  for  Soldiers’  protective  clothing,  elastomeric  adhesives  for 
composites  and  armor  design,  and  the  design  of  hypergolic  fuel  gels  to  minimize  catastrophic 
failures  of  fuel  containers  are  just  a  few  examples  in  which  viscoelastic  properties  of  materials 
are  of  critical  importance.  A  challenge  in  designing  these  types  of  systems  is  in  the  large 
parameter  space  which  needs  to  be  explored.  Many  variables  characterize  such  polymer 
systems,  including  molecular  weights,  compositions,  segmental  interactions,  and  more.  As  a 
result,  it  is  difficult  to  anticipate  the  concomitant  changes  to  material  properties  and  morphology, 
with  for  example,  temperature.  Trial-and-error  approaches  which  rely  on  laboratory 
measurements  are  costly  and  time  consuming.  Thus,  computational  material  science  has  played 
an  ever-increasing  role  in  the  design  of  polymeric  materials  (Andzelm  et  ah,  2006;  Lisal  et  ah, 
2006).  The  viscoelastic  properties  of  polymer  systems  depend  on  morphologies  and  structures 
that  span  over  length  and  time  scales  which  are  inaccessible  by  molecular  modeling.  Therefore, 
coarse-grained,  particle-based  mesoscale  models  that  retain  only  the  most  essential  features  of 
the  polymer  system  must  be  utilized. 

In  this  report,  we  outline  a  computational  protocol  that  invokes  a  particle-based  mesoscale 
method,  dissipative  particle  dynamics  (DPD)  (Hoogerbrugge  and  Koelman,  1992;  Koelman  and 
Hoogerbrugge,  1993),  to  simulate  the  viscoelastic  properties  of  various  polymer  systems  (Sen  et 
ah,  2005;  Pryamitsyn  and  Ganesan,  2006).  This  tool  will  be  invaluable  to  U.S.  Army  scientists 
in  accelerating  the  design  of  materials  with  superior  viscoelastic  properties.  In  section  2,  a  brief 
overview  of  the  dissipative  particle  dynamics  method  is  given,  followed  by  a  description  of  the 
approach  to  calculate  the  viscoelastic  properties  using  an  oscillatory  shear  technique  (Pryamitsyn 
and  Ganesan,  2006)  in  section  3,  followed  by  conclusions  in  section  4. 


2.  Dissipative  Particle  Dynamics 


The  DPD  method  is  a  mesoscale  simulation  technique  that  operates  at  time  and  length  scales 
larger  than  those  of  traditional  molecular  dynamics,  but  for  situations  that  are  inaccessible  to 
continuum  dynamics  (Hoogerbrugge  and  Koelman,  1992;  Koelman  and  Hoogerbrugge,  1993). 

A  DPD  system  is  composed  of  soft  particles,  each  representing  a  region  of  fluid,  which  moves 
continuously  in  space  and  discretely  in  time.  In  a  DPD  simulation,  the  polymer  chain  is  modeled 
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as  a  collection  of  point  particles  that  represent  lumps  of  the  chain  containing  several  segments. 
DPD  particles  are  defined  by  a  mass  mt,  r,,  and  velocity  v/,  and  interact  with  each  other  via  a 
pairwise,  two-body,  short-ranged  force  F  that  is  written  as  the  sum  of  a  conservative  force  F^, 
dissipative  force  F^,  and  random  force  F^: 

+!?+!?•  (1) 


F^  includes  a  soft  repulsion  force  F^’^  acting  between  two  particles  and  a  harmonic  spring  force 
F^'‘  acting  between  adjacent  particles  in  a  polymer  chain.  Therefore,  the  DPD  polymer  chains 
are  flexible  since  no  additional  constraints  such  as  bond  bending  or  bond  torsion  are  included. 


and  F^''^  are  given  by 
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and 
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respectively.  In  equations  2  and  3,  Uy  is  the  maximum  repulsion  between  particle  i  and  particle  j, 
=  r;  -rj,  =  |r-|,  r^,  is  the  cutoff  radius,  K  is  the  spring  constant,  and  ro  is  the  equilibrium 

spring  length.  Conservative  potential  energies  corresponding  to  F^’'  and  F^\  and  can  be 
expressed  as 


and 
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The  remaining  two  forces,  F^  and  F^,  are  given  by 


F  = 


V. 

V  y 


(4) 


(5) 


(6) 


and 


2 


(V) 
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L.  r. 

^ij  ij 
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where  (r)  and  (r)  are  weight  funetions  that  vanish  for  r  ^  r^ ,  /is  the  frietion  eoeffieient,  <j 
is  the  noise  amplitude,  =  V;  -  Vj ,  is  the  Gaussian  random  number  with  zero  mean  and  unit 

varianee  that  is  ehosen  independently  for  eaeh  pair  of  interaeting  partieles,  and  At  is  the  time 
step. 


Espanol  and  Warren  (1995)  showed  that  the  system  samples  the  eanonieal  ensemble  and  obeys 
the  fluetuation-dissipation  theorem  (in  the  limit  of  At^O)  if  the  following  relations  hold: 


®°(r)  =  [0^(r)J, 


(8) 


and 


a^=2yk3T, 


(9) 


where  T  is  the  temperature  and  ks  is  Boltzmann’s  eonstant.  (r)  and  (r)  are  typieally 
ehosen  as 


ct)°{r)-  co^[ry 


1-^ 


0 


for  r  <  /, 
for  r  >  /, 


(10) 


The  evolution  of  DPD  partieles  in  time  t  is  governed  by  Newton’s  equations  of  motion: 

^=V,(t),  (11) 

dt 

and 

(12) 

For  a  more  detailed  deseription,  see  the  original  papers  (Hoogerbrugge  and  Koelman,  1992; 
Koelman  and  Hoogerbrugge,  1993). 
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3.  Viscoelastic  Property  Calculations 


The  substantial  effort  underway  to  develop  materials  for  U.S.  Army  applications  with  improved 
thermal  and  environmental  stability  and  tunable  viscoelastic  properties  require  a  broad  range  of 
strain  rates  (10*’-10Vs)  and  elastic  moduli  (0. 1-2.6  MPa).  Therefore,  the  theoretical  methods 
which  we  propose  must  account  for  distinct  physico-chemical  processes  occurring  at  vastly 
different  strain  rates.  Both  the  measurements  and  calculations  that  attempt  to  characterize  the 
dynamic  mechanical  properties  use  a  relationship  between  the  applied  strain  and  the  resulting 
stress  in  the  material  (Ferry,  1980).  When  the  viscoelastic  material  is  subjected  to  oscillatory 
strains,  y,  of  frequency  ©,  the  stress  response,  a,  is  necessarily  cyclic  and  can  be  written  as 
a(®)  =  G"(®)  Y(®)5  where  G"(®)  =  G'  +  iwr)  and  r\  is  the  viscosity.  The  storage  modulus,  G', 
indicates  the  material’s  ability  to  store  energy  and  the  imaginary  part  of  G”(®)?  the  loss 
modulus,  ©T),  characterizes  the  amount  of  energy  lost  through  the  viscous  process.  Simulations 
can  be  used  to  determine  which  parameters  govern  the  viscoelastic  behavior  of  polymer  systems, 
delineating  the  regimes  and  frequencies  at  which  various  effects  (e.g.,  polymer  chain  length)  are 
manifested. 


The  viscoelastic  properties  are  calculated  using  a  non-equilibrium  oscillatory  shear  technique 
(Allen  and  Tildesley,  1987),  which  entails  a  simulation  with  an  additional  force  in  the  shear 
direction  along  with  time-dependent  Lees-Edwards  boundary  conditions  (Lees  and  Edwards, 
1972).  Eor  oscillatory  shear  imposed  in  the  xy-plane,  the  equation  of  motion  for  the  particle 
velocities  become 


dv.^  ^  .  d^r(t) 

^  =  Jtc+  - 


dt 


df 


(13) 


The  oscillatory  strain  can  be  taken  as  y{t)  =  A^\[- cos((»  t)] ,  where  and  ©  are  the  chosen 

values  of  the  amplitude  and  frequency,  respectively.  Eor  the  Lees-Edwards  boundary  conditions, 
the  simulation  box  and  its  images  centered  at  (x,y)  =  (±E,0),  (±2E,0),  ...  are  taken  to  be  stationary. 
Boxes  in  the  layer  above,  (x,y)  =  (0,L),  {±L,L),  {±2L,L), ...  are  moving  at  a  speed  nL  in  the 
positive  x-direction,  where  n  is  the  shear  rate.  Boxes  in  the  layer  below,  (x,y)  =  (0,-L),  {±L-L), 
{±2L-L), ...  move  at  a  speed  nL  in  the  negative  x-direction. 


Eollowing  equilibration  of  the  structure,  a  series  of  stress  trajectories  are  generated  by  imposing 
the  oscillatory  shear  conditions.  A  sample  output  for  a  diblock  copolymer  system  is  given  in 
figure  1 .  Note  that  numerical  errors  tend  to  be  higher  at  low  frequencies  requiring  averaging 
over  more  stress  trajectories  to  minimize  these  errors. 
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Figure  1.  Sample  output  for  a  model  diblock  copolymer  system  which  has  microphase 
separated  in  spherical  domains. 


4.  Conclusions 


In  this  report,  we  have  reviewed  a  computational  method  to  calculate  the  viscoelastic  properties 
of  polymeric  systems.  We  have  used  the  DPD  method  and  the  non-equilibrium  oscillatory  shear 
technique  as  implemented  in  the  code  by  Pryamitsyn  and  Ganesan  to  calculate  the  storage  and 
loss  moduli  of  the  diblock  copolymer.  Such  capability  will  be  an  invaluable  tool  for  assisting 
material  scientists  in  accelerating  the  design  of  materials  with  superior  viscoelastic  properties. 
Numerous  U.S.  Army  applications  such  as  multifunctional  materials,  coatings,  hypergolic  fuel 
gels,  and  elastomeric  adhesives  will  greatly  benefit  from  this  tool. 
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